Evidence that highly non-uniform black strings have a conical waist 
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Numerical methods have allowed the construction of vacuum non-uniform strings. For sufficient 
non- uniformity, the local geometry about the minimal horizon sphere (the "waist" ) was conjectured 
to be a cone metric. We are able to test this conjecture explicitly giving strong evidence in favour 
of it. We also show how to extend the conjecture to weakly charged strings. 



I. INTRODUCTION AND SUMMARY 

In the presence of compact dimensions massive solu- 
tions of vacuum General Relativity may take one of sev- 
eral forms including the black-hole and the black-string. 
For concreteness, consider the simplest case - d — 1 ex- 
tended space-time dimensions and a compact dimension 
of radius L, which is denoted here by the z coordinate. 
Any d > 5 may be considered, and here we specialise to 
d = 6 for numerical convenience. The transition between 
the black hole and string phases, which can be defined 
by their distinct horizon topologies, must depend on the 
sole dimensionless parameter of the problem, namely the 
ratio of the typical curvature radius of the horizon and 
L. The picture we advocate here is that black objects 
with curvature radius much smaller than L are expected 
to closely resemble a 6d Schwarzschild black hole, with 
a nearly round S 4 horizon, while very massive black ob- 
jects will have a horizon size much larger than L and will 
be extended over the z axis wrapping it completely. Such 
objects will have S 3 x S 1 horizon topology and will be re- 
ferred to as "black strings" . The simplest black string is 
a product of the 5d Schwarzschild solution with the spec- 
tator z coordinate line and we refer to it as the "uniform 
string" . Other string solutions which are z dependent 
we refer to as "non-uniform" . Gregory and Laflamme 
(GL) Hi discovered that neutral uniform black strings 
develop a perturbative z-dependent instability (tachyon) 
for masses below a critical mass m(L). Consequently the 
marginally tachyonic mode generates a branch of static 
solutions which we call the "GL non-uniform strings" . In- 
terestinglinks to thermodynamic stability were explored 
in HEIM.E1, and the instability was further studied in 

mum 

One of the authors (BK) used qualitative tools to 
analyse the phase structure of the static solutions. The 
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emerging picture consists of exactly two stable phases - 
the black hole and the (stable) uniform black string con- 
nected through the unstable GL non-uniform solutions 
in a first order phase transition. These two stable solu- 
tions have not only different horizon topologies, but the 
total topology of the Euclidean solutions is different as 
well. Thus a continuous topology change (much like the 
conifold transition in Calabi-Yau 3-folds) occurs, which 
was argued to centre around the cone over S 2 x S d ~ 3 . 

This conjecture addresses two questions. Firstly it con- 
jectures that there is only one solution for large mass ob- 
jects. This question is clearly essential for understand- 
ing how astrophysical black holes are manifested, and 
whether there are potentially observable consequences. 
Secondly, based on Gubser's higher order perturbative 
construction of the GL non-uniform strings it pre- 
dicts all strings in this branch have higher mass than the 
unstable uniform strings, and so these solutions are un- 
likely to be the end state for GL decay as suggested by 
Horowitz and Maeda I13|l33l . For another approach to 
the phase structure see |l4|. F° r implications to black 
hole uniqueness and for a description of the explosive 
transitions with Planck-scale power see 0, 0] . 

Let us recall the reason for the appearance of the cone 
and some more details. One considers a highly non- 
uniform string as well as a large black hole (a black hole 
which hardly fits in the compact dimension) and focuses 
on the vicinity of the "waist" - the minimal sphere in 
the black string, and on the "gap" between "north" and 
"south" poles in the black hole. Then one finds that 
the topology at some distance away from the waist is 
S 2 x S d ~ 3 where the second sphere comes simply from the 
spherical symmetry and the first includes the Euclidean 
time. Moreover the spatial S d ~ 3 is topologically trivial 
in the black hole (it can shrink) and correspondingly the 
S 2 for the black string. Hence a topology change anal- 
ogous to the "pyramid picture" of the conifold suggests 
itself where the cone here is over the metrically round 
S 2 x S d ~ 3 . The isometry of the cone is enhanced (rel- 
ative to the rest of the solution) and in particular the 
geometry near the waist depends on only one variable 
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- the distance p from the waist. Not only does this en- 
hanced symmetry realize the topology change in a simple 
setting, but the equations of motion actually forbid some 
more general ansatze which were tried. 

Analytic solutions are not available, despite some at- 
tempts llil Il7j and some relevant constructions in 4-d 
\wi ilyL l2Cj (where the proper size of the compact dimen- 
sion does not asymptote to a constant, and there are no 
black string solutions), which do not generalise to d > 4 
|2lj . One of the authors (TW) [H showed that in 
general the numerical method of relaxation could be ap- 
plied to gravitostatics by overcoming the conceptual is- 
sues (largely related to imposing the constraint equations 
and boundary conditions) and in particular was able to 
trace numerically the GL branch of non-uniform strings 
in the case d = 6. Large amounts of non-uniformity were 
reached (i?rnax/-R mm ~ 9) exposing the asymptotic be- 
haviour for the GL non-uniform branch including limiting 
mass and temperature. Moreover, it was found that the 
whole branch has a mass higher than the critical string 
and therefore cannot serve as an endpoint for decay of 
unstable uniform strings. The agreement with the the- 
oretical predictions of were further strengthened in 
[24| where global properties of the solutions were studied 
and shown to behave in accord with the conjecture. 

In this paper we perform a much stronger test of the 
conjecture in [Tlj by comparing the numerical geometries 
with the predicted cone. A summary of the method and 
results follows. 

First one is required to identify the transformation 
from cone coordinates to the numerical metric, and then 
to compare the two metrics. The cleanest coordinate to 
extract is p, and the "smoking gun" test involves com- 
paring the numerical Kretschmann curvature invariant 
Knum with the cone prediction, K cone , which due to 
the symmetry properties of the cone, only depends on 
p. A quantitative test shows that over a region where 
K changes by orders of magnitude agreement is found 
and K num / K cone changes only by a factor of about two, 
which we consider to be a sub- leading effect (U- Next 
we infer the azimuthal angle x horn the numerical so- 
lution. Actually, comparison of the metric gives sin(x) 
and so we must test that the maximum of the expres- 
sion for sin(x) is consistent with 1, and we see this is the 
case. Three metric functions remain to be compared, and 
indeed show good agreement. This agreement is partic- 
ularly satisfying given the low resolution of the grid over 
the relevant region, and the finite value of non-uniformity 
which can be achieved. 

In summary, we conclude that the numerical 6-d |35| 
geometry at the vicinity of the waist reveals what is es- 
sentially a cone structure, supporting the analysis of 
Sub-leading deviations from the cone remain and may 
hold further clues and refinements. We conclude with an 
argument that the waist cone geometry remains a good 
approximation if the string carries a small electric charge, 
and construct the behaviour of the electric potential near 
the cone tip. 



II. NUMERICAL STRING SOLUTIONS 

Following Gubser, we define, 

A = I(§™» -l) (l) 

where R m ax,min are the maximum and minimum sphere 
radii in the horizon. Non-uniform strings were con- 
structed non-linearly in A in 23] using the methods devel- 
oped in |22j . For technical reasons, the computation was 
performed in 6 dimensions, although it was shown that 
6-dimensional strings have the same qualitative thermo- 
dynamic behaviour as 5 dimensional strings to leading 
order in A. The numerical method uses the metric, 

ds 2 = -^- e 2A dr 2 + e 2B (dr 2 + dz 2 ) + e 2C (to + r 2 )dn 2 

(2) 

where A, B, C are functions of r, z and z = [0, L] is 
the periodic coordinate of the S 1 on which the string 
wraps. Due to reflection symmetry we actually consider 
the half period z = [0, L n ] with L n = L/2. We have Eu- 
clideanised time and now the usual Euclidean constraint 
makes r periodic, r = [0, T max ] with, 

r max - 2^ V W B ^ ) | r=0 = - (3) 

with T the horizon temperature of the string. Further- 
more A, B, C — > as r — > 00 is chosen. Choosing units 
so that to = 1, we take L — L cr u as a convenient initial 
condition for the relaxation, so L n — 2.476 and A, B, C 
decay exponentially for small A. The metric components 
are then solved by relaxation, and A ~ 4 is attainable. 
The method unfortunately is limited by resolution, and 
to proceed to higher A would require much computation 
time, and we suspect that it is more fruitful to attempt 
to improve the method itself before more detailed calcu- 
lations are performed. 

For large A the waist shrinks, compatible with van- 
ishing as A — * 00, whilst the maximal sphere radius 
Rmax appears to asymptote to a constant radius |24| . 
Other global properties of the solution were also mea- 
sured to check compatibility with the black string/hole 
transition conjecture, namely that the horizon temper- 
ature, mass, and the proper distance along the horizon 
appear to asymptote to constants for large A. Assum- 
ing the conjecture to be true, this is intuitively explained 
by the fact that at large A the geometry hardly changes 
with varying A, except near the very cone tip. As this 
region is small, it is only the very massive Kaluza-Klein 
states that participate, and these are strongly exponen- 
tially damped away from the waist, and therefore do not 
contribute much to variation of global properties. 

As in any numerical GR context, it is difficult to as- 
sess the physical significance of constraint violations. For 
small A we may compare non-linear solutions with Gub- 
ser's perturbation theory method. However, for large A 
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where we become resolution limited, it is harder to as- 
sess error. In j^] this was done by direct observation 
of the constraints, and furthermore using the first law 
to compare mass asymptotically measured with that in- 
tegrated from the horizon geometry. Good agreement 
was found, but any additional tests of numerical accu- 
racy are obviously important. The cone conjecture is a 
perfect example of such a test; at large A, if the geome- 
try is indeed a cone at the waist, we may explicitly check 
that the numerical metric does indeed reproduce this. 
That is exactly what we do here. Furthermore, if the 
conjecture is correct it may allow numerical methods to 
start relaxation about the cone geometry, for either the 
black string or black hole branch, and therefore improve 
numerical stability. 

III. TESTING THE CONE 



P 




FIG. 1: Plot of p measured from the A = 3.9 string solution. 



In 6-dimensions we may now write the Euclidean cone 
metric over the base S 2 x S 3 , 

ds 2 = dp 2 + p 2 Q (d X 2 + sin 2 X N 2 dr 2 ) + ^dQ 2 ) (4) 

where r, fl^ are the same coordinates as for the Euclidean 
string @. For the metric to be smooth away from the 
apex, x ls an angle coordinate with range \ = (0,7r). 
Furthermore, N must be chosen as 2ir/T max to ensure 
Ndr behaves correctly as an angle in the S 2 . We note, 
as emphasised in [Oj . that there are no physical free 
parameters once one has postulated this cone base. 

We may now directly compare the numerical metric 
@ and the cone Comparison of dt 2 ,d£l 2 give, 

p{r,z) = e c( "V2(l + r 2 ) 

sin X (r,z) = -;J- e A(r,z)-C(r,z) 

NV2(l + r 2 ) 

Thus p and x are completely determined and there is 
no residual coordinate freedom to complicate the com- 
parison. Knowing p and x, we may simply compute the 
remaining components of the metric Q and explicitly 
compare to (J2J. Then, 

e 2B ^ z \dr 2 + dz 2 ) = dp(r, zf + ip(r, z) 2 d X (r, z) 2 

= cl dr 2 + c2 dz 2 + 2 c3 drdz (6) 

where 

Cl = {d r pf + \p 2 (d rX f 

c2 = (d z p) 2 + \p 2 (d zX ) 2 

c3 = {d r p){d zP ) + \p 2 {d r x)(d z x)) (7) 

There are two complicating factors in executing this 
comparison. Firstly we cannot reach A — > oo but only 



A ~ 4 and thus we do not expect perfect agreement with 
the cone, as for finite A the tip will be missing. Secondly, 
numerical resolution is limited, and we only expect good 
agreement in the region near the waist. For the maximum 
A computed the r, z lattice resolution is dr = 0.025 and 
dz = 0.026. We will find good agreement with the cone is 
a region of coordinate width Ar ~ Az ^0.5, restricting 
us to a 20*20 sub-grid of the whole lattice. Considering 
large gradients are present, it is unclear how satisfactory 
this resolution is [3(|. Both of these factors add consid- 
erable uncertainty in assessing agreement. However, we 
will see that the agreement is sufficient that the cone is 
clearly a good approximation, even though it is difficult 
to make the statement more precise at this time. 



A. Comparisons of curvature invariant 

In figure ^ we plot the function p against r, z for 
the most non-uniform solution found numerically, with 
A = 3.9. We see very clearly the correct qualitative be- 
haviour of p for a cone near to the minimal sphere at 
r = 0, z — L n , where p has its minimum. The degree to 
which the cone is 'resolved' can be estimated from this 
minimum value. For the value of A plotted in the fig- 
ure, pmin — 0.27. This indicates that the cone will be 
a good approximation for p > p m i n , but p « L n , the 
characteristic scale of the compactification. 

Quantities involving only p provide the easiest numeri- 
cal comparisons. In the next section we see that comput- 
ing x from A, B, C involves further processing, and thus 
tests only involving p should be more robust. A non- 
trivial test of the cone conjecture, simply involving p is 
a comparison of the Kretschmann curvature scalar with 
the cone prediction. For the cone metric Q one finds, 

48 

K = R^ a,3 R^ a p = — (8) 
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FIG. 2: Plot showing K, the scalar curvature invariant, mea- 
sured from the string (A = 3.9) compared to the cone pre- 
diction. The ratio is also shown, and we see good agreement 
slightly away from the apex, where finite A 'resolves' the cone. 
Moving far from the cone ie. z < 2.1 or r > 0.4 the cone is no 
longer a good approximation, the actual curvature becoming 
more homogeneous in z. 

One can compute the same quantity for the numerical 
solutions directly from A, B, C and the two are plotted 
in figure [3 for the region near the waist. We also plot 
the ratio of these quantities in the same figure. We see 
that at the waist agreement is lost, as we expect because 
we are at finite A. However, moving away from the apex, 
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FIG. 3: Plots of sin \ on the z — L„ axis measured from 
string solutions with A = 1.5, 3.9. For the larger A we see the 
curve become more 'constant', and furthermore, the value is 
remarkably consistent with one, the cone prediction. 

agreement seems to be good, considering the relatively 
low numerical resolution available over such a small sub- 
region of the string solution. The ratio is close to one, 
and we also see a possible sub- leading behaviour in r, 
and less strongly in z as we move away from the 'apex'. 
Moving further from the waist we see the curvature con- 
tribution from the cone decays sufficiently that a more 
homogeneous curvature emerges, as expected for a string. 



B. Full metric comparison 

This curvature invariant indicates the size of the region 
where we expect good agreement. Now we perform a full 
metric comparison using \ as weu as P- However, use of 
X is subtle. For the cone obviously sinx = 1 for X = f/2. 
However, we are measuring sin^ directly, not Xi an d due 
to finite A resolving of the cone, and limited numerical 
resolution, we no longer can expect that sin x = 1 exactly 
at z = L n . In figure [3] we plot sin x for two values of 
A = 1.5,3.9 at z = L n , over the range of r indicated in 
the previous section to yield good agreement. Not only 
do we clearly see that for the larger A value sin x on the 
axis z = L n appears to be more constant, but moreover 
this constant is indeed consistent with one. Bearing in 
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FIG. 4: Plot of sinx computed from A = 3.9 string, including 
normalisation correction to ensure sinx = 1 at z — L n , so 
that x can De meaningfully extracted. 

mind the low resolution we regard this as remarkably 
good agreement. For the lower value of A, a constant 
behaviour is still seen at large r, again consistent with 
one, but there is a larger region near the apex where sinx 
is clearly not a constant, being approximately linear in 
r. 

In order to compute © we must compute x from the 
sinx and obviously this may lead to pathologies, as we 
see above that the sin x measured at z = L n is larger than 
one at its peak. Therefore, we normalise the sin%(r, z) 
measured by its value on the z = L n axis (which 'should' 
be one, and we have seen does well approximate one) as, 

smx(r,z) corrected = sinx(r,z)/sinx(r ,i„) 

p(r ,L n ) = p(r,z) (9) 

so we use contours of p to find the value on the axis 
z = L n to normalise by. This tidies up x at z = L n , and 
it is difficult to see how to extract x without dealing with 
this issue. 

Now we plot sin Xcorrected in figure 01 We see exactly 
the type of behaviour we would expect, with contours of 
X parameterising an angle near the apex. Then in fig- 
ure we plot the dr 2 and dz 2 components of the metric 
(cl,c2) 10 calculated using the corrected x- Accord- 
ing to the cone conjecture, these components should be 
equal to each other, and to e 2B<r,z \ which we also plot 
in this figure. We see rather remarkable consistency with 
the cone prediction. The remaining test is that the mea- 
sured drdz component (c3) should be zero, or numerically 
speaking, much smaller than the dr 2 , dz 2 components. 
This is finally plotted in figure H3 and again we see that 
the prediction is indeed true. 

IV. CHARGED STRINGS 

Uniform strings electrically charged under a Maxwell 
field can be constructed from black hole solutions of 
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FIG. 5: Plots of remaining metric components cl,c2 com- 
puted from the numerical solution with A = 3.9 using the 
'corrected' sinx- The cone predicts that cl = c2 = e 2B which 
is also plotted for comparison. Very good agreement is seen, 
given that the underlying data is on a lattice of only approx- 
imately 20 * 20 and resolution effects may give significant nu- 
merical errors. 

Einstein-Maxwell-Dilaton theory |2^, |2(| |27| . For small 
charges non-uniform deformations will exist since the GL 
instability persists. As the string horizon is an equipo- 
tential surface which is "dented" we expect the field to be 
reduced in the waist region - this is the complementary 
effect to the appearance of high fields near a "pointed 
protrusion" . Here we consider whether we can quantify 
this further, using our concrete example d = 6. 
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FIG. 6: The off-diagonal metric component c3 which the cone 
predicts to be zero. Away from the resolved apex, we see it is 
indeed much smaller than the values of cl, c2 in the previous 
plot, again being consistent with the cone ansatz. 



Unfortunately at this time we know of no generalisa- 
tion of the cone to include an electric flux. Therefore we 
simply consider putting a small charge on the cone and 
neglecting its back-reaction. Making the electric ansatz 
that the vector potential A t = (f>(p, x) with other compo- 
nents being zero, one finds, 



-d p <t> + (d 2 x 4> - cot x d x 4>) = o (10) 

The solution can be separated, 4>{p,X) — f(p)9(x) an d 
the resulting equations solved in terms of a separation 
constant fc, which we may chose to be positive. The 
radial equation is easily solved giving, 

Azk-l 



(ii) 



and since we wish the back-reaction to be small, we may 
only take solutions that are regular as p — > 0. 

Now we consider the angular equation. The horizon, 
at x — 0, 7t, must be an equipotential surface. Due to the 
separable form, we see this may either be achieved by / 
being constant (ie. k = 1), or by g — at \ — 0, tt. The 
first case yields only the trivial solution <f> = const. The 
second is more subtle, and we find the solution, 



g(x) = s\n 2 xF 



3-fc 3 + k 



2;sim 2 X 



(12) 



where F(a,b; c; x) is the hyper-geometric function. Ad- 
ditionally we add the geometric boundary condition at 
X = 7r/2, where the differential equation is singular, that 
g{ X ) t> e reflection symmetric. The hyper-geometric func- 
tion is a finite degree polynomial in sin 2 % if fc = 3 + 4n 
with integer n > 0, and hence has the required reflec- 
tion symmetry. Indeed we may Taylor expand g about 
X = giving, 



7(X) 



2^ 



'3-k s 



3+fc ' 



cosx+0(cos X 2 ) 
(13) 



which confirms that only for these values of k is the re- 
flection boundary condition obeyed. 

If we consider the exact cone geometry these modes 
with k = 3 + An clearly blow up at large p, far from the 
apex. However, for the non-uniform string, it is unclear 
what the geometry will be for large p and therefore we 
cannot specify precisely the remaining boundary condi- 
tion for the radial equation. Instead we assume that the 
leading order term is not absent in the solution, so the 
dominant contribution to the potential near the apex is 
from the lowest k, k = 3, giving, 



= g(psin X ) 2 + O(p 6 ) 



(14) 



for constant q, the higher corrections coming from n > 1. 
Therefore we see the electric field vanishes on the cone 
horizon confirming our expectation of a reduced field near 
the apex. 

To conclude, we still expect the cone to describe the 
limiting waist geometry for a weakly charged string. Ob- 
viously it would be interesting to understand how this 
picture changes as the charge is increased to extremal- 
ity where one expects the GL instability to disappear 
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